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We consider the VEV of Wilson loop operators by studying the behavior of string probes in 
solutions of Type IIB string theory generated by Nc D5 branes wrapped on an internal manifold. 
In particular, we focus on solutions to the background equations that are dual to field theories with 
a walking gauge coupling as well as for flavored systems. We present in detail our walking solution 
and emphasize various general aspects of the procedure to study Wilson loops using string duals. 
We discuss the special features that the strings show when probing the region associated with the 
walking of the field theory coupling. 
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I. INTRODUCTION 

In this paper we want to study the behavior of non-local operators of gauge theories, making use of the gauge- 
string correspondence. We are in particular interested in a specific class of supergravity solutions that are closely 
related to what goes under the name of walking in the field theory language. In the introduction we summarize the 
basic notions and ideas that will feature prominently in the paper: the treatment of Wilson loops in the gauge-string 
correspondence, the concepts of confinement and screening in gauge theories and the meaning of walking dynamics 
(with a view on its role within dynamical electro- weak symmetry breaking). 

It must be stressed that we do not know the precise nature of the field-theory dual of some of the examples we are 
going to consider in the body of the paper. The Wilson loop studied here is a very important quantity, that may help 
identifying such dual theory. 

A. AdS/CFT and Wilson loops 

According to general ideas of holography and more concretely to the Maldacena conjecture a quantum conformal 
field theory in dimension d is equivalent to a quantum theory in AdSd+i space. In general, the idea is that local 
operators in the CFT couple to fields in the AdS side, in such a way that correlators of conformal fields are related 
to amplitudes in the quantum theory in AdS, as explained in Q . For instance, since the CFT-side of the equivalence 
must contain the energy-momentum tensor T^^ among its operators, there must be a field on the AdS'-side that 
couples to it. Since this field is the graviton, then the theory on AdSd+i space must be a quantum theory of gravity. 

One can use the correspondence to study non-local operators on the CFT-side. In particular, if the field theory is a 
pure Yang- Mills theory, an example of such operator is the Wilson loop Q. These objects couple to extended objects, 
excited in the AdS side of the correspondence. The Wilson loops (path ordered exponentials of the holonomy of the 
gauge field along a curve C) are one of the most interesting observables of such a gauge theory, 

W{C) = -^TrPe' fc -^m^^" . (1) 

The loop itself and products of them provide a basis of gluonic gauge invariant operators. 

The Wilson loop along a curve C is computed in the dual string theory by calculating the action of a string bounded 
by C at the boundary of the AdS space. More concretely^, 

( W{C)) = f VFe-^^"^ (2) 

^ ' JdF{C) 

where F denotes all the fields of the string theory and dF their boundary values. A good approximation to this path 
integral is by steepest descent. The Wilson loop is then related to the area of the minimal surface bounded by the 
curve C, spanned by classical string configurations (with Nambu-Goto action Sng) that explore the bulk of AdS. 

All of this was first proposed ten years ago in In the meantime, this proposal motivated lots of developments, 
see [1] for beautiful and influential papers on this line. See also Q for a review. In particular, the ideas and techniques 
of the original AdS/CFT correspondence have been assumed to generalize to a large class of systems, and have been 
used to relate field theories that are not conformal (and hence more closely related to phenomenological applications) 
with backgrounds that are not AdS, extending the framework to what is more appropriately referred to as gauge-string 
duality. 

B. Confinement and Screening 

In its original definition the Wilson loop computes the phase factor associated to a closed trajectory for a very 
massive quark in the fundamental representation (it can also be generalized to other representations). The quark- 
antiquark (static) potential can be read from the VEV of the Wilson loop. Choosing a rectangular loop of sides 



The fundamental string is actually dual to the generalized Wilson loop, containing a term that couples the coordinates of the internal 
space with the scalars of the brane. This is due to the fact that the string ending on the brane is source of the electric field, generating 
but also of the scalars on the brane, from which is pulling. See 3] for a clear discussion of this. 
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Lqq,T, in the first approximation for large times T — ^ oo, 

(m^(C)) -e-^««^, (3) 

where Eqq is the quark-pair energy. In the hmit of large 't Hooft coupling, the steepest descent approximation 
mentioned above yields the identification 

PF(C)^ - e-^QQ"^ - e"^"^ . (4) 

The description of the Wilson loop in QCD in terms of a string partition function is not new. Well before its 
modern formulation, the ideas behind gauge/string correspondence have been used to show that the potential of a 
quark-antiquark pair separated by a distance Lqq gets a correction of the form due to quantum fluctuations of 

the Nambu-Goto action Q- 

The definition of confinement we will adopt is the following. Consider a SU{Nc) gauge theory with matter fields in 
generic representations. We decouple (make infinitely massive) all fields with non-zero N-ality, and then introduce a 
single particle-antiparticle pair of non-dynamical fields with non-zero N-ality as a test probe for the system (effectively, 
the pair dynamics is quenched). We then compute the work needed to separate the particle-antiparticle test pair up 
to a distance Lqq. If the work approaches Eqq ~ ctLqq for large separations, the theory is confining^. The quantity 
CT is a representation dependent constant, the string tension. 

C. Walking technicolor 

Walking technicolor [8J is a framework within which the phenomenological difficulties of dynamical electro-weak 
symmetry breaking might find a very natural and elegant dynamical solution, thanks to the fact that, in contrast 
to QCD-like theories, the guidelines provided by naive dimensional analysis are violated. This is because of large 
anomalous dimensions controlling the dynamics over a large regime of energies. The word walking refers to the fact 
that the new dynamics is strongly coupled over a large range of energies, where its fundamental coupling exhibits a 
/3-function which is anomalously small in respect to the coupling itself. A behavior of this type is expected in theories 
which flow onto strongly-coupled IR fixed-points, and it is reasonable to assume that it persists also when such IR 
fixed points are only approximate, though in this case a degree of ambiguity as to the meaning of approximate invites 
some caution. 

Besides being affected by the usual calculability limitations due to the strong coupling (as in QCD and in QCD- 
like technicolor), the walking dynamics itself makes this framework very hard to work with. New, non-perturbative 
instrument are needed in order to understand the (effective) field theory properties of a walking theory. Very recent 
years saw a lot of progress towards a better understanding of walking dynamics both from the lattice Q and thanks 
to ideas borrowed from the gauge-string duality. See for instance [Ifll for a list of references focused on the precision 
electro- weak parameters. 

In this paper, we will use the word walking to refer to backgrounds in which there exists a interval in the p 
radial direction over which the geometry is determined by some coupling that shows an evolution with the radial 
coordinate that is anomalously slow. This agrees with the standard definition of walking, where the beta-function of 
the fundamental gauge-coupling is, over some range of energy, anomalously small in comparison with what expected 
from the strength of the actual coupling itself. Not all couplings show this behavior: this also agrees with standard 
definitions, in the sense that relevant operators must be present in order to drive the RG fiow away from a possible IR 
fixed point, and scale invariance is present only in the sense that in the walking region the relevant couplings are so 
small that their effect can be neglected. However, our definition is significantly less restrictive that the usual one: we 
do not require the presence of an actual fixed point in the flow, and correspondingly we do not have an approximately 
AdS background, nor can we recover an AdS background by dialing some parameter to some special value. 

D. General idea: string-theory as a laboratory for walking dynamics. 

Motivated by the difficulties described in the previous subsection, in ref. a more general program is proposed, 
based upon gauge/string correspondence in order to go beyond the low-energy effective field theory description. The 



^ Another equivalent way of defining confinement is by computing the VEV of the Polyakov loop, that if vanishing indicates a confining 
theory. Also, the perimeter law of a 't Hooft loop indicates confinement. 
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proposal is to study the dynamics of theories that yield walking behavior, but that are not necessarily related to the 
electro-weak symmetry. In short, one would like to use string theory as a laboratory in which to study the general 
properties of walking dynamics by itself, in isolation from its complicated realization within an explicit model of 
dynamical electro- weak symmetry breaking. In ref. [ll|, it is shown that in the context of Type IIB string theory on 
a background generated by a stack of Nc D5 branes, there exists a very large class of solutions to the background 
equations for which a suitably defined gauge coupling exhibits the basic properties of a putative walking theory. The 
running of the gauge coupling flattens over a large range of intermediate energies, but restarts at low energies, until 
the space ends into a (good) singularity in the deep IR, so that no exact IR fixed point exists. The fact that this 
is not a walking technicolor theory (there is no electro-weak symmetry in the set-up, and hence no mass generation 
in the usual sense), together with the large- A^c expansion, yields the advantage that we avoid the complications due 
to mixing of weakly-coupled and strongly-coupled properties of the theory. For instance, the spectrum of the spin-0 
sector of the theory can be studied, and has been studied [l2j, yielding remarkable surprises. 

In this paper we take a further step in this direction, by studying the behavior of the Wilson loop in backgrounds of 
this class. As we will explain, we can use the techniques developed in the context of gauge-string duality, by studying 
the background with a probe string. In particular, we will study Wilson loops in the dual walking QFT. For technical 
reasons that will be explained in the body of the paper, in order for this program to be carried out we will also need 
to generalize further the class of backgrounds in These new solutions have been already introduced in [l^. We 
explain here in deeper detail how to generate them, characterize them, and relate them to the literature. 



The paper is organized as follows: we set up notation and introduce a set of important ideas in section ITTl revising 
the bibliography and adding important new ingredients and derivations. Then we apply these to well-established 
examples of gauge-string duality in section IIIIl providing a simple and compact set of exercises that are intended to 
yield some guidance in the following sections, in which the dynamics is far from well-understood. Section ITVl presents 
our new walking solution and a study of the dynamics of the Wilson loop as a function of the length of the walking 
region. Section FVl studies the results derived in section HIl when applied to background that encode the dynamics of 
fundamental fields. We then conclude in section IVll 



In this section we present general results for Wilson loops, computed using the ideas of Some of the results here 
have been derived long ago (see for example [l3|), but our approach will be different and some new and useful points 
will be specially emphasized. 

We study the action for a string in a background of the generic form 



We assume that the functions {gtt, gxx, gpp) depend only on the radial coordinate p. By contrast, gij for the internal 
(typically compact) space can also depend on other coordinates. Whatever are the internal coordinates, they will play 
no role in what follows. This is because we will choose a configuration for a probe string that is not excited on the 0* 
directions, hence in what follows, we will ignore the internal space ^. 



E. Outline 



II. GENERAL THEORY 



gttdt'^ + gxxdx^ + gppdp^ + gijdO^dO' . 



(5) 



A. Equations of motion 



The configuration we choose is. 



t = T, x = .T(cr), p = p{(t). 



(6) 



and compute the Nambu-Goto action 




(7) 



^ Strings or other objects that extend in the internal space filling part of it can be treated as an effective string, analogous to the one we 
are studying. If these objects are allowed to vibrate in the internal space, then a generalization of the present treatment should be done. 
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The induced metric on the 2-d world-volume is Gq^ = g^,ydaX^dpX^ ^ where 

Grr — ~gtt, Ga 



''da' 



Defining for convenience f{p)^ = gttgxx, 9ip)^ — 9ttgpp, the Nambu-Goto action is 



T 
2na' 



2tt 



T 
2na' 



2ir 



daL . 



(8) 



(9) 



Notice that we consider the situation in which the string does not couple to a NS B-field. 

We first compute the Euler-Lagrange equations from Eq. ([7]) and then we specify them for the ansatz in Eq. ([5]). 
We get that for the {t,x,p) coordinates the eqs. of motion read, respectively, 



dr 
da 



-1 
.1 

ri 

I 
1 



{fx'^+g'p'^) 



0, 



fx' 



0, 
1 

I 



(10) 



{x^ff + p'^ga'), 



where X' = '^'^J^J^ for any function X . 

The first equation in (jlOp is solved because we assume a background metric independent of time (we consider the 
system at the equilibrium) . The second equation in (jlOp is solved if the quantity inside brackets is a constant (that 
we call C), which implies that. 



dp -^f 
da \da / 



_(J2 



(11) 



One can check that the third equation in (|T0| is solved if Eq. (fTTj) is imposed. Then, we need to work with just one 
equation. Defining 



we write it as 



dp , dp 

^ = ±^v;,,(p)o- = ±K//(p). 



(12) 



(13) 



Another way to arrive to the last version of Eq. ([T^ is to consider a restricted ansatz for the string configuration 
[i = r, X = a, p = p{a)] and use the conserved Hamiltonian derived from Eq. ^ to get an expression for p[a) that 
is precisely Eq. ([T^ . 

The kind of solution we are interested in can be depicted as follows: a string that hangs from infinite radial position 
at a; = and drops down towards smaller p as x increases. Once it arrives at the smallest p compatible with the 
solution, namely po, it starts growing in the radial direction up to infinite p where x = Lqq, see Fig. ([l])^. This means 
that in the two distinct regions x < Lqq/2 and x > Lqq/ 2 the equations of motion will differ only in a sign 



X < 



Lc 



X > 



Lc 



We can now formally integrate the equations of motion 

x{p) 



TOO dp 

Jp K„M 



X < 



(14) 



(15) 



For future reference we refer to the lowest end of the radial coordinate as po ■ 
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: on 



FIG. 1: Setting of the string. 



or more compactly 



dr 



Veffir) 



(16) 



In what follows we will use only one of the solutions in Eq. (jl4[) unless explicitly noted. 



1. Boundary conditions 



We need to specify the boundary conditions for the string in Eq. This is an open string, vibrating in the bulk 
of a closed string background. Following the ideas of \^ , we add a D-brane at a very large radial distance where the 
open string will end. This string will then satisfy a Dirichlet boundary condition at p — >■ cx). This means that for 
large values of the radial coordinate ^ must vanish. The only way of satisfying the equation of motion in Eq. (fTTj) 
for p — ^ oo, given that the left hand side has to be non vanishing, is to have a divergent Veffip)'- 



lim Veff{p) 
p^oo 



(17) 



This implies that there are restrictions on the asymptotic behavior of the background functions [f{p),g{p)\ in order 
for the string proposed in Eq. © to exist. We will come back to this in the following sections. Before proceeding, a 
brief digression is needed. When studying Eqs. (jl3p - (|17p the reader may find unsatisfactory that the restriction we 
have proposed above, namely 



dp 

dx 



Veff\p^^^CO 



(18) 



p— >oo 



looks dependent of our choice of the radial coordinate. This is actually not the case, because we could rewrite this 
restriction in a more covariant form in the following way. We define a couple of vectors ^ in the x, p directions. 



da 



da''"' 



(19) 



^ We thank Johannes Schmude for the discussions that lead to this. 
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compute the ratio ^ of their norms and impose that this is divergent on the surface V on which the string has to 
satisfy the Dirichlet condition: 



yxx y ^„ ) 



V 



oo. 



(20) 



Combining this with Eq. (|13p. and evaluating at the boundary 



9pp 

9 XX 



fj2 



(21) 



This last expression is free of coordinate ambiguities, as it comes from operating with invariants (norms of vectors). 
It is however easier to work within a specific choice of coordinates, which we will do in the body of the paper. As we 
will see explicitly, this occurs in the examples we will study below. 



2. Turning points 

Once Eq. (jl7|) is satisfied the string will move to smaller values of the radial coordinate down to a turning point 
Po where the quantity ^(po) = Oj i- e. 14// = 0. In principle, there could be points where Ve// vanishes because of 
either isolated zeros of /(p), or diverging point of g{p). However, we will not consider this kind of inversion points, 
since we are interested in solutions of the equations of motion that allow the string to probe the entire radial direction. 
Hence the turning point can be placed in any possible po: with po < po < oo (where po is the end of the space). Thus 
we will restrict ourselves to forms of T4// where the inversion point is given by imposing C — f{po). Furthermore, in 
the next section we will also impose that the envelop of the string is convex in a neighbourhood of po, hence insuring 
that gauge theory quantities like the separation between the pair of quarks and its energy will be continuous functions 
of Po. It is then clear from Eq. (jl3p that Ve//(p) controls not only the boundary condition at infinity, but also the 
possibility for the string to turn around and come back to the brane at infinity. 



B. Energy and Separation of the QQ pair. 

We now follow the standard treatment for Wilson loops summarized in [gJ . If our probe-string hangs from infinity, 
turns around at a point po as described above and goes back to the D-brane at infinity, we can then compute gauge 
theory quantities, like the separation between the two ends of the string, which can be thought of as the separation 
between a quark-antiquark pair living on the i'-brane and coupled to the end-points of the string. And we can 
compute the Energy of the pair of quarks, that we associate with the length of the string (computed along its path 
in the bulk). Both of these quantities will be functions of the turning point po. 

The standard expressions that we will use can be derived easily. Indeed, for the QQ separation we only need to 
compute / dx. To calculate the energy of the QQ pair, we compute the action of the string and substract the action 
of two 'rods' that would fall from infinity to the end of the space ^. The results are @, 



Lqq{po) = 2/(po) 



' 9{z) dz 

Po 



■oo 



EQQipo) = f{po)LQQipo) + 2 / ^VPiz) - P{Po)dz - 2 / g{z)dz . (22) 



Po J \ I ■J Po 



As discussed above, the constant C defined around Eq. pT|) must be taken to be C = /(po)- Using Eq. ([T^ we can 
rewrite 



°° dz 



Lqq{po) = 2 / — ^ (23) 

/po ^eftiz) 



Notice that these 'rods' need not be strings, it is just a way of rcnormalizing the infinite mass of the quarlcs. 
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As discussed in section III A 11 we must impose that for large values of the radial coordinate T4// diverges. This, 
however, is not enough to ensure that the integral above receives a finite contribution from the upper end of the 
integral, we then have to require that for large radial coordinate Vef / diverges at least as 



p— >oo 



(24) 



with /3 > 1. Then, the quantity Lqq can be finite or infinite depending on the IR (p ^ po = po) asymptotics of the 
background, meaning that we consider the turning point at the end of the space at the lower end of the integral. If, 
expanding around the turning point po, we have 



P->Po 



(25) 



then the separation of the quark pair is infinite for 7 > 1 and finite for 7 < 1. We will see examples of both behaviors 
in the following sections. 

Another reasonable condition that we could impose is that po is the first zero of Veff for which the string has 
positive convexity (a minimum). This can be easily expressed as a condition on the background functions. In fact, let 
us assume that near po the effective potential in Eq. (|T2|) behaves as Veff = k{p — poV ■ Then, in a neighbourhood of 
Po, we have 



dp 

dx 



±Veff = ± {k{p - poT) + 0{p - por+' 

Veffip)- 



dx"^ dx'' dx' dp dx <-JJ\f^ 

by iteration and induction we obtain (we choose only one of the branches of Eq. (|14p ). 

^ = (^^//(p)d^)> = «"n^=o [j(7 - 1) + 1] (P - Po)(^-^)"+^ + o[p poY''-'^ 



-l)ri+2 



(26) 



(27) 



In order to have a minimum we impose that the first non vanishing derivative at p = po is an even derivative (its 
value has to be positive). So, for even n we need (this result will not depend on the choice of branches above), 



7 = 1- 



1 



(28) 



from this it follows that 7 cannot be bigger than one. This means that the only possibility of obtaining a string that 
can be stretched up to infinite distance will come from the case 7=1 (contrast this with the result of [i3|)- The 
solutions with infinite length will have all even derivatives vanishing at pp. Had we considered non- integer corrections 
to the leading term in Eq. (|26p . so that near po the function 



^eff = Kl(p- po)^ + K2(P- Po) 



7+e 



(29) 



for very small values of e we would have had a formula like the one in Eq. (|27p where the exponent is the same, but 
the coefficient changes. For large values of e, we have, to leading order, the same result as in Eq. (I?7l) . The reasoning 
given above applies. 

These constraints ensure that there exists a unique trajectory for the string as function of po. It should be also 
noticed that the profile of the string is not analytic function and in particular in po it turns out to be a function. 



1. Some exact results. 



Next, let us derive an expression of the energy Eqq as function of the inter-quark separation Lqq. For this purpose, 
it is useful to introduce the function K[x] = ,J, ^ . The expression for Lqq{pq) in Eq. ([2^ reads^ 



Lqq(po) = lim 2 r 



/(po) 



dz . 



(30) 



Some of these expressions have been derived in past collaborations with Angel Parades. 
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Performing the derivative 



dLc 



dpo 



J[Po) 



lim 2 



dnj< 



/(^) 



Now, we use the following identity, 



9poK 



fiPo) 



fiPo) 



f{z)fiPo) 
f'{z).f{po) 



fiPo) 



integrating by parts and after some algebra, we get 



or 



dpo 

dLqq 
dpo 



' lim 2dp„\ogf{po) 

pi— >-oo 



29polog/(po) lim 



gjpi 

f'ipi 



-K 



fiPi) 



pi 



dz K 



fiPo) 



pi 



Po 



Po 



/(Po) 



dz K 



fiPo) 



lipi) 



Similarly, we compute the derivative for the energy function finding that 



dpo 



fiPo 



dLc 



dpo 



and we get 



dE, 



QQ 



dL. 



QQ 



fiPo). 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



Notice that (after integration) Eq. (j36p is an exact expression for Eqq in terms of Lqq , the non triviality residing 
in the function po{Lqq), this expression was also derived in [13] ■ Also, Eq. (pG)) yields the force between the quark 
pair. It should be stressed that we are assuming that we will be able to invert po as function of Lqq. Whenever this 
cannot be done, the solution would be valid only in the regions of monotonicity of Lqq(po)- 

Let us comment briefly about the existence of cusps in our strings. If near the lower point of the string (typically 
situated in the end of the geometry) the quantity Veff{p) diverges, then there will be a cusp in the shape of the 
string. This can be easily understood using Eq. and an expansion of Veff{p) of the form discussed in Eq. p^ . 
If Veffip) cc {p — pa)'^ with 7 < 0, one can integrate the equation 



P - Po oc 



(1-7) 



'QQ 



(37) 



which implies that the shape of the string p{x) is not analytic at — f^- 



C. Leading and subleading behaviors. Inversion points 

Let us focus on the case that, near the end of the space (in the IR), Vgff ~ (p — Po) (the string stretches up to 
infinite length Lqq{po) -^00). The subleading corrections to the Energy of the quark pair [l3| read 

Eqq = f{po)Lqq + K + 0(e-l-l^«« (log Lqq)") . (38) 

This formula can be obtained as an expansion around po ~^ PO: or equivalently expanding Eq. (j36|) around Lqq — > 00. 
Notice that no power-law corrections or Luscher-terms appear in these corrections: they would appear if 7 > 1 in 
Eq. (|28|) . something that we ruled out. Contrast this with the result of [l^. Because Eq. ([28]) and the discussion 
around it rely only on the generic properties of Ve//, such power-law corrections can only descend from higher- 
derivative corrections to Eq. (O, which cannot be repackaged into the form of Eq. p^ . 

In the following sections, we will discuss these subleading behaviors in various examples and precisely find the 
coefficients that apply to each particular case. Before proceeding, a few comments are due in order to avoid ambiguities. 

The partition function of the string is studied by using two different expansions, in a' and gg. First, the Nambu- 
Goto action is characterized by the string tension, and one can think of expanding any physical quantity (correlation 
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function) in powers of a'. This procedure is effectively equivalent to quantizing the (particle- like) excitations of the 
string, and in this sense ol corrections can be associated with loops of the dynamics of the string modes. One can also 
rephrase the results in terms of a classical effective theory, in which the oi corrections are encoded in higher derivative 
corrections to the Nambu-Goto classical action. By doing so, one can associate the resulting Eqq to the expectation 
value of a Wilson loop in the dual theory, and hence interpret the generalization of Eq. (l38l) in terms of the actual 
quantum corrections to the quark-antiquark potential at large distance, for a confining string, as done for example 
in [l5| . As explained above, this procedure yields power-law corrections to the leading-order result Eqq cx Lqq. This 
is noi what we are doing in this paper. All the results we obtain are at the leading-order in a', the calculations being 
completely classical and based on the Nambu-Goto action. In particular, this explains why we do not find a Luscher 
term. 

We are going to truncate at^ the leading-order also the expansion in g^^ , which controls processes where the string 
breaks. This is also done in |l5j . where the quantum corrections that are computed are the ones in a', but the whole 
analysis treats the string itself as effectively free. Not only quantum corrections related to the Qs expansion, but the 
whole many-body nature of string-theory is neglected in this way. This is a very important limitation: strings that can 
stretch up to infinite separation for the quark pair can be treated in this way, while theories in which fragmentation 
and hadronization take place via string-breaking are accessible to our approach only up to a scale smaller than the 
scale of breaking itself (real world QCD, if it admits a string dual description, should fall into this class in which 
corrections must be included and actually dominate the dynamics from the scale of breaking and below) . 

We will see in the following examples in which Eq. (I^5|) holds with 7 < 1. This will yield a finite value for the 
maximal displacement in the Minkowski directions between the quark-antiquark pair. Studying the large-Lgg limit 
would necessarily require the inclusion both of quantum effects in the o! expansion, but also the effect of string 
breaking, which might be decoupled in the large quark-mass limit. 

Another important comment is the following: if the function Lqq{pq) is not monotonic, then it is not invertible, and 
Pq{Lqq) is multivalued. If this is the case, it will happen that for some given Lqq we can find that also Eqq{Lqq) 
is multivalued. Among all the possible values of Eqq for a given Lqq (all of which satisfy the equations of motion) 
we will refer to the lowest one as the stable solution, while the others will represent either excited or unstable states. 
The Van der Waals liquid-gas system provides a nice realization of these excited and unstable states, and we report 
its treatment in Appendix [B] 

The existence of inversion points (or extrema of Lqq{pq)) implies that the first derivative vanishes. Using 

the expression in Eq. ([M)) . 



rPi 

lim i dz K 



pi- 



/(Po) 



OA^tK)- . HO. (39) 



/'W; f'iPl)Veffipi) 



If it happens that f'(pfyv^ff(pi) ~^ ^ large values of pi, then a necessary condition for the existence of turning 
points is that the integral in Eq. (p4)) vanishes, or more simply that the integrand changes sign. Notice, nevertheless, 
that this is a not a criterium of practical application except for some particular easy backgrounds as we will see in 
section |Vl 

To close this section let us stress that given the constraints we have specified along this section, the shape of the 
string that we are proposing as solution is the only allowed shape, since there can be only one minimun of p{x). 

III. SOME WELL-UNDERSTOOD EXAMPLES 

In this section we illustrate the points of the previous section in a set of very well known examples. All of them are 
string-theory backgrounds that are simple enough that the dynamics of the probe string can be treated analytically. 
We will write the 10-dimensional background in string frame. 

A. The case of AdSs x 5"^ 

This is certainly the main example, where the conjecture was originally proposed [l'| and the ideas for Wilson loops 
described in the introduction first developed Q. After a rescaling of the radial coordinate, the metric reads, 

p „2 j-)2 -, 

ds^ = a' ^dxl 3 + —dp^ + R^dnl , (40) 
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where = ^AirggNc is dimensionless, p has dimensions of inverse-length, and the constant a' in front of the metric 
ensures that when we compute the functions f,g defined below Eq. (jS]) and plug into Eq. ^ the factors of a' will 
cancel (this will happen in the examples discussed below also). Hence we have, 



fip)' = ^"^ yeff^^^P'-pt (41) 



pIR' 



One can check that the functions of the background respect all the constraints we imposed in section [IT] regarding the 
boundary conditions and convexity. We can exactly integrate Lqq{pq) 

Lqq(Po) = 2/ dp ^ (2^)i-ii^ (42) 

Jpo PWP -Pt Po^ijV 

Since it is possible to invert the relation we can then write pq{Lqq), and using Eq. p6p 

dLQQ LqqL (j)* LqqL ( j)* 

in agreement with the result of [3]. Of course, the separation for the quark pair diverges for po (the end of the 
space). Also notice that in this case, the expression of Lqq{pq) is invertible. There are no corrections to Eq. 

A few words of comment might be useful to a reader who is not familiar with this set-up. The Yang-Mills coupling 
of the dual N — A field theory is defined as gy^ = '^'^9s, so that the (dimensionless) curvature R'^ is proportional to 
the 't Hooft coupling in the CFT. The supergravity approximation holds when i?** ^ 1. 



B. Witten-Sakai-Sugimoto Model 



This model, based on DA branes wrapped on a circle with SUSY breaking periodicity conditions ^TB\ received a 
great deal of attention thanks to the observation by Sakai and Sugimoto [17] that the introduction of Nf ^ Nc flavor 
D8 branes as probes allows to construct a model where a peculiar realization of chiral symmetry is spontaneously 
broken. The metric reads 



d^ 



3/2 



dxi^ 



F{p)dxl 



R 



3/2 



dp^ 

Hp) 



p''d^^ 



(44) 



where is the coordinate on the circle, = Trg^Va' and F{p) 



Notice that in this case the gauge 



coupling gyM = {'^'^Y9sVa' is dimensionful, so that p has dimensions of inverse-length, as in the previous subsection, 
and so does A. The relevant functions are, 



[C^^f{po?^a''§ 



12 p" 

12 p^ 

p^-Ai" 



' (p3_A3)(p3_^3) 

[p^Rf 



(45) 



hence. 



Lqq = 2{poR) 



3/2 



dp 



v/(p3-pj5)(p3_A3) 



(46) 



One can explicitly perform the integral above, but the result, in terms of special functions, is not very illuminating. 
To get a better handle on the underlying dynamics, we consider the case in which po = A. In this case, we get 



Lqq{po) = - 



Gpl 



12 arctan 



'2p + Pa 
VSpo 



log 



PPo 



pI 



{p-poY 



(47) 



that we can see diverges logarithmically for p = pQ. So, the string here again, has infinite length as is expected in a 
dual of a confining field theory. If po > A then the separation of the pair is computed from the integral of Eq. (j46p 
and turns out to be finite. 
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On the other hand, it is possible to iteratively invert the relation between Lqq and po for K and hence for 

Lqq -> oo 



Po(iQQ) = A + 4V3e A + ... 

and from Eq. (j36p we get for the Energy of the pair Eqq in terms of the separation Lqq , 

^QQ(iQQ), = ^QQ p -^(^ A) + ... 



(48) 



(49) 



C. D5 branes on S 



In this example the system consists of Nc D5 branes that wrap a two cycle inside the resolved conifold preserving 
four supercharges (A/" = 1 in four dimensions) [1^. After a geometrical transition takes place, a background dual to 
this field theory contains a metric, dilaton and RR three form. In this case gyM = (27r)^5sa'- We quote only the part 
of the metric relevant to this computation (after a rescaling of the Minkowski coordinates is done) 



dx\ 3 + dp^ + dsj 



The relevant functions read, 



(50) 



l{p) = f{p)^e^gsa'N, 



e^'^Q sinh(2p) 



(51) 



2, /pcoth(2p)- 



The integral defining the separation of the quark pair cannot be evaluated explicitly, but we can check that the upper 
limit of the integral gives a finite contribution (it goes as J°° p^e~Pdp), while from the lower extremun of the string 
reaches the end of the space {po — >■ 0) we get 



dp 



pa 



=0 ye2^opM~ Po^O 



lim log(/9o) — > oo. 



(52) 



indicating that (like in the Witten-Sakai-Sugimoto example discussed above ) strings that reach the end of the space 
correspond to an infinite separation between the quark pair, which in turn indicates the absence of screening (expected 
from the QFT that contains only fields with zero N-ality). 

We compute the subleading terms in the Energy of the pair in terms of the separation to obtain, 



EQQ^e^^'^^LQQ + Oie-'^^'^^) 
as above a clear sign of a confining dual QFT. 



(53) 



D. Klebanov-Strassler Model 



Certainly, the Klebanov-Strassler model is the cleanest example for a dual to a four-dimensional field theory that 
confines in the IR and approaches a conformal point in the UV (modulo important subtleties) fl9l |. Here again, we 
will quote only the relevant part of the metric. 



ds' = h{p)-'/'dxl, + Mp)1/2,4/3^p!^ ^ ^^2 



6if(p)2 



with 



^3^^^ ^ sinh2^^ Mp) = 2^/^(^^.a'^c)^^-«/^ 
2smh (p) 



; coth X — 1 



sinh 



2 (sinh2x- 2a;)^/^. 



(54) 



(55) 
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The functions /, Ve// read, 



In this example again, using the asymptotics for the various functions, it can be checked that when the string 
approaches the end of the space po = the separation of the quark pair diverges fogarithniically, as in the two 
previous examples. Also, the energy of the pair reads 

Eqq = f{0)LQQ + ©(e-^^^Q) (57) 
with flQ ~ (g a'jv^e-'t/^)^ ^ 1.1398... As above, it is clear that the model shows confining behavior. 



4/3 _ VW^)K{p) 



IV. WALKING SOLUTIONS IN THE D5 SYSTEM, UNFLAVORED. 

This section is devoted to the specific case of a class of solutions to the D5 system that exhibit walking behavior in 
the IR, in the sense that a suitably defined gauge coupling becomes almost constant in a finite range of energies [ll| . 
We remind the reader about the set-up, based on the geometry produced by stacking on top of each other D5- 
branes that wrap on a S*^ inside a CY3-fold and then taking the strongly-coupled limit of the gauge theory on this 
stack (that leaves us in the supergravity approximation for this system) . We introduce the class of solutions of interest 
and apply the formalism of the previous sections to these solutions. As we will see, these solutions do confine, in the 
sense defined earlier on, but also show a remarkable behavior in the walking energy region, leading to a phenomenology 
resemblant a first order phase transition. As a result, the leading-order, long-distance behavior of the quark-antiquark 
potential is linear, but the coefficient is very different from what computed in section fill CI 



A. General set-up. 



We start by recalling the basic definitions that yield the general class of backgrounds obtained from the D5 system, 
which includes [18] as a very special case. We start from the action of type-IIB truncated to include only gravity, 
dilaton and the RR 3-form F: 



1 



Sub = -p^r- I d. xy/-g 



We define the SU{2) left-invariant one forms as, 

LUi — cos "tpdO + sin sin 9d(fi , W2 = ~ sin -tpdO + cos "0 sin 9dip , uj^ = dip + cos 9d(fi . 



(58) 



(59) 



and write an ansatz for the solution 23] assuming that the functions appearing in the background depend only the 
radial coordinate p, but not on x nor the 5 angles 0, 9, (p, </>, "0 (in string frame): 



a 9s 



e29(p) g-^Hp) 
— — ((wi + a{p)d9f + (w2 - a[p)su\9d'ff) + - '~- ' — ^-'■-^^ 



-(cja + cos 9d(py 



- (a)i + b{p)d9) A (0)2 - b{p) sm9dLp) A (^3 -I- cos 9dip) + 



b'dpA {-d9 A cji + sm9dip A UJ2) + (1 - Kp)'^) s.m9d9 Adip Auj^ 



(60) 



The system of BPS equations can be rearranged in a convenient form, by rewriting the functions of the background 
in terms of a set of functions P{p), Q{p), Y{p),T{p), cr^p) as [2l| 



4e 



2/1 



P cosh T — 



= Pcoshr- Q, 



2k 



AY, a 



P sinh T 

P cosh T — 



Neb = a. 



(61) 
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Using these new variables, one can manipulate the BPS equations to obtain a single decoupled second order equation 
for P{p), while all other functions are simply obtained from P{p) as follows: 

Qip) = (<3o + Nc) coshr + Nc{2pcoshT - 1), 

sinhT(/?) = . ^ — — — , coshT(p) = coth(2p - 2po), 
smh(2p — 2po) 

P' 

Y{p) - ^, 



e*"^" cosh(2po)2 



(P2 _Q2)ysinh2r' 

CT = tanhr(Q + TVc) = — r-r-T;^ „ - x ■ 62) 

smh(2p — 2po) 

The second order equation mentioned above reads, 

In the following we will fix the integration constant Qo = ~-^ci so that no singularity appears in the function Q{p)- 
We also choose po = for notational convenience, together with a'gg = 1. For our purposes, it is also convenient to 
fix Se^"^" = 1. With all of this, the functions we need for the probe string are: 



.2/ N 20 / sinh^(2p) 

f (P) - = \JjJSTZrQ2^^ (64) 

g'ip) = \p'fip), (65) 



B. Walking solutions. 

The simplest and better understood solution of the system described in the previous subsection is defined by 

P = 2N,p. (67) 

It belongs to class I, and it has already been presented in section IIII C[ where is shown the study of the behaviour 
of the Wilson loop^. 

For this background it is possible to show that 

fiO) = (68) 

^ ' NcV2Na 

and that for = /^(0)i expanding around p ^ 

In particular, on the basis of what we already discussed around Eq. psp this means that Lqq and Eqq diverge for 
C2 ^ /2(o). 

In [11], it was observed that there exists a class of well-behaved solutions for which a suitably defined (four- 
dimensional) gauge coupling exhibits a walking regime, meaning that for a long range in the radial coordinate p/^ < 



see Appendix |X] for the definition of class I and II and for more details about the solutions to the equation for P. 
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p < the running becomes very slow, and the gauge couphng effectively is constant. These solutions depend on two 
free parameters c and a, and can be obtained recursively, by assuming c large and expanding in powers of iVc/c, so 
that 



P{p) = ^ci-"Pi_„. (70) 

n=0 

with 

1/3 

Pi(p) = (cos^ a + sin^ a(sinh(4p) - Ap)) . (71) 

In particular, for c/Nc — ^ +oo, the solution is very well approximated by P ~ cPi. 

The solutions of the form in Eq. ([70|) belong to class II in the language introduced in Appendix [A] This type 
of solution is not suited for the present study, because of the exponential behavior of P and P' at large-p, which 
is not compatible with the boundary conditions for the string in the UV as required in Eq. (ITTl) . (Equivalently, 
the presence of high-dimensional operators dominating the dynamics in the far UV renders the study of the probes 
problematic. Analogous problems arise when studying the spectrum of excitations of the background 12*] and/or of 
probe fields [13] •) However, we are mostly interested in what happens in the IR. We are hence going to construct and 
study a different class of solutions, which can be thought of as a generalization of Eq. (|70l) . with UV-asymptotics in 
class I. Such solutions can be expanded for small-p, yielding, 

P{p) = Co + kzCQp + k^cop" H ^ -^^^ h • • • , (72) 

with Co and k^ the two integration constants. Notice how this expansion does not contain a term linear in p. This 
means that in the co limit one does not trivially recover Eq. (I67p . 

In order to build numerically the solution, we start by expanding Eq. (|63p . by assuming that the solution can be 
written as 

P{p) = Pip) + efip) , (73) 

and replacing in Eq. 

^ Fo + eFi + 0{e^) . (74) 
Because P is an exact solution, Fq — 0. Hence one finds a new equation, that now is linear: 



(/(p)(- cosh( 4p) + 4psinh(4p) + 1) - 2psinh ^(2p)/'(p)) 
- 4 sinh(4p)p + cosh(4p) - 1 
8((l-4p)/(p)+p/'(p)) 



Ap-l 



rip) (76) 



In the last step, we approximated the equation by assuming that p ^ 0. The resulting equation can be solved exactly, 
yielding, in terms of hypergeometric functions: 

/(p) = e-4''y4^ciL/(^^,|6p-0+e-4py4;^C2L(-^,i 6p-0 . (77) 

Asymptotically, neglecting power-law corrections, this means 

/(p) ~ Cie-4P + c^e'P , (78) 

implying that consistency of the perturbative expansion in Eq. (TM)) enforces the choice C2 — 0. Indeed, there are 
no asymptotic (in the UV) solutions that behave as e^^ . Allowing for C2 ^ would imply that the solution is not a 
deformation of P, but rather that the expansion in Eq. ([75)1 breaks down, and the solution (if regular) falls in class 
II. The procedure of allowing for a small component ci 7^ can be interpreted as the insertion of a small, relevant 
deformation, which does not affect the UV-asymptotics, but has very important physical effects in the IR. 

We have now obtained an important result: at least up to large values of p, there exists a class of solutions that 
approach asymptotically the P solution. We cannot prove that such solutions are well behaved all the way to p 0. 
However, we can use this result in setting up the boundary conditions (at large-p) and numerically solve Eq. 
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towards the IR. By inspection, these sohitions are precisely the ones we were looking for. They start deviating 
significantly from P below some > 0, below which P is approximately constant. We plot in Fig. [2] two such 
solutions, with ~ 4 and ~ 9, together with the P solution for the same value of Nc. We also plot in Fig. [3] the 
functions appearing in metric (e^^, e^^, e^*^, 0) for the same solutions. Notice the behavior of e^^ for p — > 0, but also 
the fact that the dilaton (p is finite for p — >■ 0. 

A short digression is due at this point. The reader should be aware that in the presented case we are likely working 
with a singular background, as the function e^^, a warp factor inside the internal space, is divergent at p = 0. 
The problem of singularities in gauge-string duality has a long history. Surely it is better to work with non singular 
backgrounds, but it is possible to obtain interesting information also from singular manifolds. The procedure, devel- 
oped along the years, to decide whether one is a dealing with a "good singularity" or not, consists in analyzing all 
the physical observables and showing that no track of the singularity can be found in them. One specific suggestion 
is to look at the behavior of the gu component of the metric (in the present case, the dilaton) and make sure that 
it is finite. Many are in literature the examples of such a class of backgrounds. To corroborate the idea that our 
background falls in this class, we decide to investigate not only the value of the Wilson loop but also two invariants of 
the metric, namely the Ricci scalar and a suitable contraction of the Ricci tensor. As can be seen in Fig. |4]no trace of 
singularities can be found in these two observables for our background, in support of the idea that the backgrounds 
we are considering are indeed acceptable. 




2 4 6 8 p 10 



FIG. 2: The numerical solutions for P{p)/Nc used in the analysis as an example, for A'^c ~ 100. The three solutions correspond 
to the P case with p, — 0, and to two new numerical sohitions with, respectively, p, ~ 4 and p, ~ 9. The numerical solutions 
can be plotted up to p ~ 150, but in the following we will truncate them at pi — 30. 



C. Probes: numerical study. 



We set-up the configuration of the string by assuming that its extremes are attached to the brane at p = pi ^ 1, 
and treat this as a UV cut-off. In the numerical study and in the resulting plots, we used pi — 30. The string is 
stretched in the Minkowski direction x — x{p), with 2^(pi) — for convenience. We vary the integration constant 
> /^(O). For each choice of C we define po as V^yy (po) = 0. In this way, the coordinates of the string are {x{p), p), 
where 

dr 

The Minkowski distance between the the end-points of the string is hence Lqq — 2x{pq). For the energy, because in 
the numerical study we do not remove the UV cut-off, rather that Eq. ([22]) . we use the unsubtracted action (setting 



17 




FIG. 3: The functions {e^^ , e^^, e^*, 4>) appearing in the metric for the same solutions as in Fig. [2] computed rescaling P — >■ P/Nc, 
and Q -> Q/N^. 




FIG. 4: The (10-dimensional) Ricci scalar R and the scalar = RmnRpqq'^'^ g'^^ , plotted as a function of the radial 
coordinate p, for several numerical solutions in the class discussed in the body of the paper. Each curve corresponds to a 
different value of p*. Notice that both scalars are finite in the p — >■ limit. 



Tl{2-Ka') — 1) evaluated up to the cut-off: 




= 2/ (80) 

The numerical results obtained for the three different solutions P are shown in Fig. [S] and Fig. [51 
Let us first focus our attention on the P case of Eq. (|67)) . The deeper the string probes the radial coordinate 
(smaller values of po), the longer the separation Lqq between the end-points on the UV brane, in agreement with our 
criteria and results of sections |TT] and IIIIl 

Two regimes can be identified: as long as po > PiBi then Lqq varies very little with poj while for small po, further 
reductions of po imply much longer Lqq. The scale pm ~ C'(l) is the scale in which the function Q changes from 
linear to approximately quadratic in p, and is also the scale below which the gaugino condensate is appearing (the 
function h{p) in the background is non-zero). This result is better visible in the upper panel of Fig. [5] The dependence 
of Lqq on Po is monotonic, but shows two very different behaviors for po < piR and po > pir, respectively. The 
transition between the two is completely smooth. 
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FIG. 5: Upper panel, the radial coordinate po of the middle point of the string as a function of Lqq. Middle and lower panel, 
the energy Eqq as a function of the quark-antiquark separation Eqq{Lqq). The three solutions in Fig. [2] are used, with the 
same color-coding. 
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The physical meaning of this behavior is well illustrated by studying the total energy Eqq of the classical configu- 
rations, as a function of Lqq and of po (see the middle and lower panel of Fig. [5]). One sees that for small Lqq, the 
energy grows very fast with Lqq, until a critical value beyond which the dependence becomes linear. We have already 
studied analytically this behavior, which can be interpreted in terms of the linear behavior of the quark-antiquark 
potential obtained from the Wilson loop in agreement with the discussion of sections HIl and Hill The energy is also a 
monotonic function of pQ. 

By comparing with the solutions that walk in the IR, one sees that a very different behavior appears. Starting 
from the upper panel in Fig. [SJ one sees that as long as po > p*, the dependence of Lqq from po reproduces the P 
case. Beginning from such large po, we start pulling the string down to smaller values of po, and follow the classical 
evolution. Provided we do this adiabatically, we can describe the motion of the string as the set of classical equilibrium 
solutions we obtained in the previous section. Going to smaller pq, Lqq increses, and nothing special happens until 
the tip of the string touches po — P*- At this point Lqq = Lmax ■ From here on, the string can keep probing smaller 
values of po only at the price of becoming shorter in the Minkowsi direction (smaller Lqq). Another change happens 
when Lqq = Lmm, at which point the tip of the string entered the bottom section of the space, po < piR. From 
here on, further reducing po requires larger values of Lqq. Asymptotically for po — >■ 0, the separation between the 
end-points of the string is diverging, Lqq — oo. 

Even more interesting is the behavior of the energy (lower panel of Fig. [Sj: for very short Lqq, and again for very 
large-^QQ, it is just a monotonic function, but for a range Lmm < Lqq < Lmax there are three different configurations 
allowed by the classical equations for the string we are studying ^. One of the three solutions (smoothly connected 
to the small-^QQ configurations) is just the Coulombic potential already seen with P. The highest energy one is 
an unstable configuration, with much higher energy. The third solution (smoothly connected to the unique solution 
with Lqq > Lmax) rcproduccs the linear potential typical of confinement. Notice (from the lower panel of Fig. [5]) 
that the solution at large-Lgg is linear, but has a slope much larger than what seen in the P case. This co-existence 
of several disjoint classical solutions is expected in systems leading to phase transitions (see Appendix [B]). The 
instability /metastability /instability of the solutions can be illustrated by comparing Fig. [5] with Fig.[TU]in Appendix IB] 
This is just an analogy, and one should not push it too far. However, identifying the pressure P, volume V and Gibbs 
free energy G as Lqq -(r^ P, pq ^ V and Eqq o G, one sees obvious similarities. In particular, there is a critical 
distance Lmm < Lc < Lmax at which the minimum of Eqq is not differentiable. 

In order to better understand and characterize the solutions we find, it is useful to look more in details at the 
shape of the string configurations, focusing in particular on the middle panel in Fig. ^ in which we plot the string 
configuration that solves the equations of motion for various values of po, on the background with p^ ~ 4. Consider 
those strings that penetrate below p*. Besides having a shorter Lqq, and higher Eqq than those for which po > P*, 
this strings show another interesting feature. They start developing a non trivial structure around their middle point, 
that becomes progressively more curved the further the string falls at small p. Ultimately, this degenerates into 
a cusp-like configuration, which disappears once pq approaches the end of the space. Notice that, as a result, the 
three different solutions for Lmin < Lqq < Lmax have three very different geometric configurations. One (stable 
or metastable) configuration is completely featureless, and practically indistinguishable from the solutions in the 
background generated by P. The second (stable or metastable) configuration shows a funnel-like structure below 
p*, and then the string lies very close to the end of the space. The third (unstable) solution presents a highly 
curved configuration around its middle point All of this seems to be consistent with the fact that in the region 
PiR < p < p* the background has a higher curvature, and as a result the classical configurations prefer to lie either 
in the far-UV or deep-IR. It must be noted here that the Ricci curvature is indeed bigger in this region, but that it 
converges to a constant for p — > 0. 

D. Comments on this section. 

We conclude this section by summarizing here three important lessons we learned. 

First, we are able to derive the linear Eqq{Lqq) expected from confinement. This linear behavior emerges for 
Lqq ^ Lc, Lc being the critical distance between the end-points of the string, of the order of the distance Lqq 



Notice that the qualitative behavior of the solutions with p« ~ 4 and p* ~ 9 are identical. However, we were not able to follow 
numerically the solution with larger p*, and hence the plots show only two such solution. The existence of the third is assured by the 
fact that this background must yield a confining potential. 

This should not be interpreted literally as a cusp. The classical solutions we found arc always continuous and differentiable, and this 
structure disappears once the string approaches the end of the space. 
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FIG. 6: The strings in (a;, p)-plane, obtained with various choices of > /^(O). Top to bottom, the three numerical solutions 
corresponding to increasing values of p* 
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computed for the string the tip of which reaches p/jj. This fact provides a physical meaning for the scale pir, which 
clearly shows in the (scheme-dependent) background functions and in the gauge coupling 

In particular, we can conclude that the walking solution we looked at is dual to a confining theory, but a very 
different one from the one of the non-walking solution P. This is signaled by the fact that the leading order behavior 
of the walking solution is related to the value of Vg/ / in po > while the non- walking is related to the zero of Vg/ / 
in po }t 0) ^-nd hence the presence of this two different solutions has to be related to the existence of two different 
scales in the background. 

A limitation of this formalism is that, as we explained at length, it does not allow to study the sub-leading corrections 
to the linear behavior, which would require a treatment in which a'-suppressed quantum corrections are included. It 
would be very interesting to calculate these corrections, which might provide useful information as to the nature of 
the dual theory. 

The second important lesson we learn is related to the second scale p* appearing in solutions with walking behavior. 
Again, this scale appears in P and as a consequence in all the functions in the background, including the (scheme- 
dependent) gauge coupling. The main point is that, provided p* > piR, p* has a very important physical meaning: it 
separates the small-Lgg regime, where the background and all physical quantities are the same as in the original P 
solution, from the large-Lgg regime, in which the dual theory is completely different from the dual to the P solution. 
The scale p* is not just a scheme-dependent fluke effect: its value is somehow related to the coefficient of the linear 
leading behavior of the quark-antiquark potential. 

The third lesson we learn is that in the region piR < p < p^, the dynamics being more strongly coupled than 
elsewhere, classical solutions are unstable. The string configurations that are stable are those that either do not reach 
p*, or those that go through this region only in order to reach the region near the end of the space, where the string 
can lie up to indefinitely large separations (in the Minkowski directions). This is a very interesting result, that might 
be related with what found in [l2| , where it is shown that for large values of p^, the spectrum of scalar glueballs on 
these backgrounds splits into a set of towers of heavy states and some light state, separated by a large gap. This is 
going to be studied elsewhere [23|. 

V. WILSON LOOP IN A FIELD THEORY WITH FLAVORS 

In this section we consider the effects of fundamental (flavor) degrees of freedom on the Wilson loop, using back- 
grounds that encode the dynamics of fields charged under a flavor group. As discussed in the Introduction, we expect 
that confinement does not take place, instead the theory will screen. We will observe the existence of a maximal 
length, requiring the string to break. As we explained, we are not including effects hence the pair-creation will not 
be accessible to the description we are giving here. In [13] these effects are taken into account by constructing the 
screened solution explicitly. 

The construction of string backgrounds where the effects of flavors is included was considered in a large variety of 
models, see [l^. We will concentrate on the models developed in [2^. We follow the treatment and notation of plj . 
as done in the previous sections. Other authors have studied effects similar to the ones described in this section by 
using different string models, see [27|-[3l|. 

At this point, let us stress that, what we will do in the following is to use as probe a closed string in a classical 
theory with no gs correction. Such an observable cannot describe the physics of the broken QQ pair in the QFT 
sector; It will have a good overlap with the state describing the QQ pair only for small separation distance among the 
quarks. If the distance among the quark is increased another configuration, that will mimic the breaking of the pair, 
will become more energetically favorable. On the other hand if we decide to measure only the Wilson loop we will 
see no signal of such a breaking. Conversely if a cusps appears in the profile of the string used to measure the wilson 
loop, this signal has to be read as the loss of validity of the approximation imposed on our theory. While this signal 
cannot be interpreted as corresponding physical quantity (the distance at which the cusp appears has nothing to do 
with the scale of breaking of the QQ pair), the observation can give a suggestion about the range of validity of our 
approximation. The results of this section should be understood as the application of the formalism of section |ll] to 
the case of flavored background, with the proviso of the validity of the Nambu-Goto approximation and the comment 
made above about the screening phenomena. Aside form this, the numerical solution Fig. [7] below is new material 
included here. 

In order to flnd the background solutions, we will need to solve a differential equation that is a generalization of 



A change of scheme in the string picture corresponds to a redefinition of the radial coordinate. 
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the one discussed in the previous section, Eq. (|63| : 

-P' + Q' + 2Nf P'-Q' + 2Nf 



P" + (P' + Nf] 



P-Q 



P + Q 



^ - 4coth(2p) 



0, 



Q{p)= \ ^ (2pcoth(2p)-l). 



(81) 



The relation to the functions that appear exphcitly in the background is given in Eqs. (3.18) and (3.19) of |2l| . 

There are various known solutions to Eq. ((8T|) . We focus on the so called type-I solutions that are known only as a 
series expansion. For small values of the radial coordinate (p — >■ 0), the expansion is given in Eq. (4.24) of ^ while 
for large values of the radial coordinate {p -> oo) is given in Eqs. (4.9)-(4.11) of (2]| : 



P{p^Q)^P„-N}p+-c'P^p' 

Nf Nf{Nf~2Nc) 
P(,^oo) = Q + iV.(l + ^+ 



(82) 



In the Nf — )■ limit this can be matched with the expansion Eq. ((72|) . Below, we find it useful to have the IR {p 0) 
asymptotics for the functions 



p40 



e40o r 47V. 107V2 



Po 



and in the UV (p — )■ oo) 



p2 
^0 



„40 



o4p 



P 



6^'= ^ 2c3 



e^'^ ^ 1. 



P2p' - 2PoNfp' 



We plot in Fig. [7] a set of numerical solutions, which reproduce the asymptotic behaviors in Eq. 

25 r 

Pip) 




(83) 

(84) 
for several 



FIG. 7: Numerical solutions of P(p) for various values of Nf/Nc- 

different values of Nf/Nc- The existence of these solutions smoothly joining the asymptotic behaviors has been 
assumed in [2l|. Here we are explicitly showing that this assumption is correct. Numerically, we learn that (provided 
Po is not too small) the function P is convex, as shown in Fig. [T] Let us move to the study of string probes in these 
backgrounds. 
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In the light of the discussion of section III Al we form the combinations 



f2 _ 

J — QttQxx 



-.2^M- V8e^^"sinh(2p) ^2 ^ ^2,(po) 



/ — yttypp — c — 



eff 



k{p) 



(85) 



We will choose for convenience the integration constant 00 such that 8e'*"^° = c^Pqi hence we will have (l){p = 0) = 0. 

Following the discussion below Eq. ((23)) and using eq. (|84p . the contribution to Lqq of the integral for large values 
of the radial coordinate goes like 



Lqq{po ^- 00) - / 

which is finite. Similarly, using Eq. (j83p . the lower-end of the integral contributes with 



dp 



Lqq{po 0) 



dpy/p, 



(86) 



(87) 



Po^O 



which is itself finite. In this case we can see that in Eq. (|25|) the quantity 7 = — 5, according to the discussion around 
Eq. ([^5)) implies a finite Lqq . Contrast this with the examples studied in section IIIII and [Vj 

The fact that near p = the quantity Vg// ^ p^^ implies, using Eq. (|13p . that near the IR there must be a cusp-like 
behavior. A plot of the shape of the string makes this clear, see Fig. [8] where we string probing a background with 
Nf/Nc = 1.2. In contrast with the example studied in section ITVl these backgrounds have a divergent Ricci scalar at 
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FIG. 8: Shapes of the string for various choice of po 



p = 0, in agreement with the fact that only for p^ — Q the string presents a cusp. This is the only example presented 
in this paper with a true cusp in the string profile. 
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A. The case Nf = 2N^. 

Finally, we turn our attention to the very peculiar, limiting case with Nf — 2Nc. The solution is described in 
Eq. (4.22) of the first paper in It corresponds to 

Q = ^N,^^-^^,P^N, + ^nTTQ^, (88) 

and the functions f{p),g{p) read in this case 

f =9^ = e2'^"+2^ (89) 

where now the radial coordinate is defined in the interval —oo<p<oo. The explicit computation of Eqs. (|22p . now 
for the lower end of the integral pq — > — oo, gives 

Lqq = ^g,a'N,7:, Eqq = (90) 

We illustrate in Fig. [S] the behavior of the string on this background. Notice how the length of the string is fixed 
Lqq — TT, in units where a'gsNc = 1, irrespectively of how deep into the radial direction the string probes the 
background. Another interesting fact is that the quantity of Eq. p9p vanishes identically in this background, meaning 
that (ILqq / dpo — for any po, implying that it is always possible to find a solution with any given po E^nd with the 
same Lqq, as Fig. [9] clearly indicates. The field theory dual to this background is quite peculiar. On the one hand, the 
numerology Nf = 2Nc brings to mind the situation in which the theory becomes conformal. Nevertheless, we should 
not forget that this is a background constructed with wrapped branes, so, clearly a scale is introduced (the scale set 
by the inverse size of the cycle wrapped). It may be the case that once the singularity is resolved one finds in the deep 
IR an AdS^ space. But the theory should remind that in the UV conformality is broken, hence developing the scale 
A^'^ = a'gsNc, that is typical of six dimensional theories. Notice that the only allowed separation, according to the 
calculations above is precisely tt in units of this scale and the energy is zero. These results do not need to be generic 
for all solutions with Nf — 2Nc, but for this particular one discussed above. Note also that in this case, even when 
the background is singular in the IR, the string does not become cuspy, as the case studied in the previous section. 
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FIG. 9: The string hanging on the special background in Eq. (|88() . 



VI. SUMMARY AND CONCLUSIONS 

Let us summarize what we learnt in this paper and emphasize the reasons that motivated our study. In section 
|TT]we recovered many well-known results about holographic Wilson loops, but at the same time introduced a certain 
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amount of formalism and some new results that allow us to make a systematic study of many qualitative features of 
string probes dual to Wilson loops in field theory. For example, we introduced the quantity V^fj in Eq. ([T^ which 
dictates the behavior of the string near the boundary of the space (p — >■ cx)) and tells us if the Dirichlet boundary 
condition is attained. At the same time, Veff near the IR tells us about the possibility (or not) of stretching the 
probe string indefinitely. Assuming certain characteristic power-law behavior of the function Veff near the IR, we 
were able to derive a condition on the power ensuring the presence of a minimum (or turning point). In that same 
section, we derived an exact expression of the relation between the energy and the separation of the quark pair. We 
also found an integral expression giving a necessary condition for the existence of turning points in the string probing 
the bulk and the condition for the existence of cusps. 

In section Um we applied the formalism described above to some well-understood examples and made some comments 
about the absence of Luscher terms in the Eqq — Lqq relation. 

The material in section IIVI was derived with the following motivation: in the paper a new background was 
proposed to be dual to a QFT with a walking regime (for a particular coupling). Nevertheless, aside from an argument 
based on symmetries given in that paper, it could have been the case that this walking region is just a fluke of the 
coordinates choice that dissapears under a simple diffeomorphism (in the dual field theory, gauge couplings and beta 
functions are scheme dependent). To clarify if there is any real physical effect in the walking region we computed 
Wilson loops in the (putative) walking field theory. In order to do this, we needed to find a new walking solution that 
allows the string to satisfy the boundary conditions discussed in sectionUT] The new solution was found numerically as 
a perturbation of the solution in 18] . We gave an interpretation of this new background in the dual field theory as the 
appearance of a VEV for a quasi-marginal operator. When studying the dynamics of strings in this new background 
we observed that there are various physical effects associated with the presence and length of the walking regime: the 
value of the Ricci scalar, the separation between the quark pair, the relation between the energy and the separation of 
the pair, etc. In conclusion, the walking regime has physically observable effects, hence it can not erased by a change 
of coordinates (conversely, it is not an effect of a choice of scheme in the field theory). Recent experience seems to 
suggest that whenever we deal with a system with two independent scales, the phenomenology of the Wilson loop will 
be similar to what we described in section HVl 

Finally we moved to the study of the Wilson loops in field theories with fiavors. Here we benefitted from some work 
done in the past, where the asymptotic behavior of the background functions was given. We constructed numerically 



the full solution in terms of a convenient formalization of the problem described in |2l| . We applied the material of 
section [n] to this case and checked that the probe string was behaving as expected (screening). 

We believe that this paper clarifies a considerable number of new and interesting points. Certainly, it pushes 
forward the idea of using string theory methods to study models of walking dynamics. This may become important at 
the moment of modelling the mechanism electroweak symmetry breaking. On the phenomenological side we provide 
a set of tools that can be applied to other backgrounds and may be even used in bottom-up approaches to the dual 
of QCD. 

It would be nice to find new models of walking dynamics and apply the ideas in this paper to study and compare 
features. It would also be interesting to apply the formalism developed here to some of the examples of backgrounds 
dual to field theories with flavors, for which certain aspects of Wilson loops were studied [13] - [Slj . 
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A: UV asymptotic solutions. 



In this appendix we present a short digression about the asymptotic behavior of the solutions to Eq. (j63p , in order 
to make contact with a possible field theory interpretation of the results that have been presented in the section HVl 
Generically, one can expand all the functions appearing in the background, and associate the integration constants 
appearing in this expansion with the insertion in the UV of the dual field theory of operators with scaling behavior 
determined by the p-dependence of the corresponding term in the expansion. This is reminiscent of what is done in 
the case of backgrounds that are asymptotically AdS, though in the present case the procedure is less solid, and the 
results should be taken with caution. 

We can start from the function Q. 

Q ^ 7V,(2p-l + 0{e-^f)) , (Al) 

where we dropped factors as powers of p in the exponential corrections. We will do so in rest of this analysis, the 
logic of which will not be affected. 

If one thinks that the two behaviors correspond to the insertion of an operator of dimension d and of its 6 — d 
dimensional coupling in the (six-dimensional) theory living on the I?5 branes, this means that this expression should 
be written as 

g - z'' + z^-'^ , (A2) 

with z proportional to a length scale. By comparison with the expansion of Q, one is lead to the identification 
p = — llogz. This result agrees with what was done by comparing an appropriately defined beta- function to the 
NSVZ beta function in ,22||. One can interpret the two terms in this expansion as the deformation of the theory by 
the insertion of a marginally relevant operator of dimension d ~ 6 — e, with coupling of dimension Q — d ^ e. Notice 
that the coefficient of the 0{e~*P) correction is very small, and has a sizable effect only at very small values of p. This 
is not due to the fact that we chose a particular value of the integration constant Qq, in order to avoid the arising 
of a pathology in the IR. Allowing for Qq would yield the same expansion, and the modification of the coefficients 
would be such that only for p < piR would the sub-leading correction have an effect. The scale pm ^0(1) has an 
important role in the present study. 

Another important quantity in the background is 

smh(2p) 

which by means of the previous identification is equivalent to 6 ~ z'^. This can be interpreted as the VEV of a 
dimension-3 operator (customarily identified with the gaugino condensate in the dimensionally-reduced 4-dimensional 
low-energy theory). This is a relevant deformation. Because this is a more relevant operator than the VEV mentioned 
above, it always dominates in the IR over the sub-leading correction to Q, which can be ignored. Notice also that 
b{p) is practically vanishing for p > pjR. This provides an intuitive explanation for this scale: it is the scale at 
which the gaugino condensate emerges dynamically. It also suggests that the expansion of Q, rather than identifying 
an independent operator, might be interpreted in terms of the insertion of an operator loosely corresponding to the 
square of the gaugino condensate. 

There exist two different classes of UV-asymptotic solutions for P [2lj : 

P - 2Ncp + 0{e-^P) (class I), (A4) 

P - 0{e^/^P) + 0{e-^/^P) + 0{e-^/^P) (class II). (A5) 

In class I, there is only one integration constant, in the 0{e^'^'') term. Notice that P has the same leading and sub- 
leading components as Q. However, the sub-leading correction depends on a free parameter: the corresponding VEV 
can be enhanced in such a way that there be a range of p over which the dynamics is dominated by this deformation, 
over the deformation present in b{p), while the latter will become important only at very small values of p. 

Solutions in the class II are rather different. The independent coefficients appear in the ©(e'*/^'') and 0{e~^/^P) 
terms, and the former cannot be dialed to zero independently of the latter (see [2l| for details). Using the same 
identification between p and z, the leading order component of P scales as z~^, and can be interpreted as the 
insertion of a dimension-8 operator in the six-dimensional theory. It is somewhat natural to think that it is related to 
a gauge coupling is six-dimensions. The presence of sub-leading corrections that scale as z^, and z* suggest that the 
gravity field P should not be interpreted as a simple operator in the underlying dual dynamics, and that some caution 
should be used. One might want to interpret the z"* as the insertion of the VEV of a four dimensional operator. But 
it could as well arise from the combination of the coupling scaling as z~^ and the same VEV of marginal operator 
scaling as z^ that is present in class I. All of this should not be taken too literally, but provides some guidance in 
what we do in the body of the paper, when constructing the class of walking solutions we are interested in. 



27 



B: Van der Waals gas. 



Here we summarize some aspect of first order phase transitions that plays an important conceptual role in the body 
of the paper. By way of example, we remind the reader about the classical treatment of the van der Waals gas, in 
terms of its pressure P, temperature T and volume V N moles of particles, by means of the equation of state 



P 



NRT N^a 



V ~hN 



where R, b, a are constants. 

In Fig. [TUl we plot one isotherm. The condition for stability of the equilibrium {d^F/dV^^, 
satisfied in some region. This implies that a phase transition is taking place. 



(Bl) 



(If )^ > is not 





FIG. 10: The pressure P as a function of the volume V (left panel) and the Gibbs free energy G as a function of the pressure 
P(right panel) for the same isotherm curve. 



In order to understand what the physical trajectory followed by the system at equilibrium is, we plot in Fig. llOl the 
Gibbs free energy G = G{T,P) = G{P) for the same isotherm. From this plot, one sees that the system evolves on 
the path ABCOQ, where C = A^, in such a way that for every choice of P the Gibbs free energy G is at its minimum. 
The evolution is smooth along ABC (gas phase: \dP/dV\ is small), but at C = iV the free energy is not differentiable, 
signaling that a first-order phase transition is taking place. In the (P, V) plane the system runs along the horizontal 
line (constant P) joining C and N . This explains the Maxwell rule introducing a curve of constant pressure that 
separates two regions of equal areas above and below it, delimited by the original isotherm. Afterwards, the evolution 
follows smoothly the curve NOQ (liquid phase: \dP/dV\ is large). 

While the trajectory ABGNOQ follows the stable equilibrium configurations, it is possible to have the system 
evolving along the path GDE or LMN . Both these paths represent metastable configurations, because {dP/dV)rp < 0. 
Indeed, these metastable states can be realized in laboratory experiments. For example, in a bubble chamber or in 
supercooled water, the metastability is exploited as a detector device, because small perturbations induced by passing- 
by charged particles are sufficient to drive the system out of the state and into the stable minimum. The evolution 
along the path EFHIL is completely unstable and not realized physically (it is a local maximum, as clear from the 
right panel of Fig.[lOl and by the fact that {dP/dV)j. > 0). 

The analogy with the examples in the main body of the paper is apparent. Notice for instance that the pressure P 
in this system as a function of the volume V behaves as a non monotonic function, hence there are inversion points 
in the curve (the points P, L in figure [T0|) . in the same sense in which we discuss the presence of inversion points in 
the body of the paper. 
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